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ABSTRACT 

The existence of submillimeter-selected galaxies (SMGs) at redshifts z > 4 has recently been con- 
firmed. Using simultaneously all the available data from UV to radio we have modeled the spectral en- 
ergy distributions of the six known spcctroscopically confirmed SMGs at z > 4. We find that their star 
formation rates (average ~ 2500 M Q yr -1 ), stellar (~ 3.6xl0 11 M Q ) and dust (~ 6.7xlO 8 Af ) masses, 
extinction (Ay ~ 2.2 mag), and gas-to-dust ratios (~ 60) are within the ranges for 1.7 < z < 3.6 
SMGs. Our analysis suggests that infrared-to-radio luminosity ratios of SMGs do not change up to 
rcdshift ~ 5 and are lower by a factor of ~ 2.1 than the value corresponding to the local IR-radio 
correlation. However, we also find dissimilarities between z > 4 and lower-redshift SMGs. Those 
at z > 4 tend to be among the most star-forming, least massive and hottest (~ 60 K) SMGs and 
exhibit the highest fraction of stellar mass formed in the ongoing starburst (~ 45%). This indicates 
that at z > 4 we see earlier stages of evolution of submillimeter-bright galaxies. Using the derived 
properties for z > 4 SMGs we investigate the origin of dust at epochs less than 1.5 Gyr after the 
big bang. This is significant to our understanding of the evolution of the early universe. For three 
z > 4 SMGs, asymptotic giant branch stars could be the dominant dust producers. However, for 
the remaining three only supernovae (SNe) are efficient and fast enough to be responsible for dust 
production, though requiring a very high dust yield per SN (0.15-0.65 Mq). The required dust yields 
are lower if a top-heavy initial mass function or significant dust growth in the interstellar medium 
is assumed. We estimate lower limits of the contribution of SMGs to the cosmic star formation and 
stellar mass densities at z ~ 4-5 to be ~ 4% and ~ 1%, respectively. 

Subject headings: galaxies: active - galaxies: evolution - galaxies: high-redshift - galaxies: ISM - 
galaxies: starburst - submillimeter 
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1. INTRODUCTION 

Submillimeter-selected galaxies (SMGs) are among 
the most powerful starburst galaxies in the Universe. 
Most of them have been found at redshifts 1.5-3 
(IChapman et al. [20 05). Their complex selection crite- 
ria ( Blain et al.l 2004?) . in particular the requirement of 
a radio detection to obtain a precise localization, make 
it difficult to discover the very high redshift tail of 
SMGs. This was addre ssed using deep, high resolution 
observations of SMGs Jfono et al. 2006; Tacconi et al.l 



2001 [200l fWang et al.ll2007l 120091: lYounger et alJl2007L 
2008ahbl . l2009bl lal: iDannerbauer et al.l 120081: ICowie et al l 
20091 ) The existence of SMGs at z > 4 has re- 




At these redshifts the age of the Universe is < 1.5 Gyr, 
which enforces the need for careful analysis of the 
timescales for formation of stars and dust. The impor- 
tant question is if supernovae (SNe), or asymptotic gi- 
ant branch (AGB) stars, or some other sources are re- 
sponsible for production of dust residing in these galax- 
ies. Locally, dust is pred ominantly f ormed by evolved, 
post-main-sequence stars (|Gehrd[l989l), but the si t uation 
may be different at high redshifts. iDwek et al.l (|2007l ) 
claimed that only SNe can produce dust on timescales 
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less th an 1 Gyr, but it has been shown bv lValiante et all 
(2009) that AGB stars begin to dominate dust produc- 
tion over SNe as early as 150-500 Myr afte r the onset 
of star formation (see also lSloan et al.ll2009f ). SN-origin 
dust has been claime d to be present in z ~ 6.2 quasar 
(jMaiolino et al.ll2004f ) and in two gamma -ray burst host 
galaxies at z ~ 6.3 (jStratta et al. 120071 but se e other 
interp retation of their da ta in lZafar et al.ll2010| ) and at 
z ~ 5 (jPerlev et al.ll2009D . 

In order to understand the formation of SMGs and 
their evolution through cosmic time, it is also important 
to compare high- and low-redshift SMG samples. This 
may help to constrain when their stars were formed. 

SMGs at z > 4 are also suitable to study the in- 
frared (IR) - radio correlat ion. This remarkably tight 
corre lation, found locally (Hcl ou et al.1 1985; C ondon 
1992), was studied a t redsh ifts 
20021 iGruppioni et all 120031: 



< 3.5 dGarrett 



Kovacs et al 
20071: 



>pioni et 
.1 120061: 



Bovle et al 
Vlahakis et all [2007; 



Apple ton et all 12004 



12007E 
lYang et al.l 



Marleau et alJ 
2007 



Beswicket all l2008t lHainlinel 120081: llbar etail 1200 



Sauna et alTT20 08: Ga m et al J l2009t iMich alowsk i et all 



iMurphv et all 120091: ._ 
Seymour et alJ 120091:" 



Murphvl 120091: iRieke et al I 
.Younger et alJ ~l2009b; 
Sargent et al.l 120101) . No significant evolution of the 
correlation was found up to these redshifts, but 
SMGs seem to form a correlation by their own offset 
toward higher radio luminosities dKoyacs et al.l 120061: 
Michalowski et al.l 120101 : IMurphv et all 120091 : IMurphv! 
2009). The only sign of evolution was reported by 



Ivison et a l. (2010) based on stacking analysis of 24 /im 
selected galaxies, though possibly interpreted 
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TABLE 1 

Photometry Detections of z > 4 SMGs 







^obs 


Flux 


Error 




SMG 


2 


w 


(MJy) 


(A*Jy) 


Ref. 


LESSJ033229.4-275619 


4.760 


0.597 


0.072 


0.006 


Coppin09 


LESSJ033229.4-275619 


4.760 


0.652 


0.215 


0.020 


Coppin09 


LESSJ033229.4-275619 


4.760 


0.771 


0.350 


0.013 


Coppin09 


LESSJ033229.4-275619 


4.760 


0.905 


0.405 


0.015 


Coppin09 


LESSJ033229.4-275619 


4.760 


3.600 


2.860 


0.060 


Coppin09 



Note. - The measured fluxes (Column 4) are given 
with the corresponding la errors (Column 5). The data for 
SMMJ163555. 5+661300 have been corrected for lensing magnifi- 
cation of a factor of 5.5. (This table is available in its entirety in 
a machine-readable form in the online journal. A portion is shown 
here for guidance regarding its form and content.) 

Ref e rences. - ICapak et al.l 12004. 120081): IGiavalisco et al.l 
i200lh; Wang et al.l 120041. I2007i, l2009l)TTFrayer et al.l 120061) 



Iono et al.1 120061); IPope et all 1 2006); Danncrb auer et al.1 (2008) 
Greve et al. ,' (2008); P erera et al . (2008); Schin nereT et all (EiuM 
Casey et all 120091) ; lc"hapin et al.l 120091) ; ICoppin et all 120091) 
Daddi et ail (2009bl la()- TDevlin et all 120091 ) ; IKnudsen eHd] 120101 ) 



TABLE 2 

Photometry Upper Limits of z > 4 SMGs 



ICapak et aT] (120041, |2008|); IGiavalisco et~aT 
12004 120071. 120091); IFro 

) 

120091) 1 

IDevlin etaTJ 12009); 



(2004); WangetaU 12009. 120071. I2009D ; IFraver et al.l 
Iono_et_^Q (2006); Poge^t^lJ^ ( 2006 ) ; Danncrbauer et al 
Greve et al. (2008); Perera et al. (2008 
Casey et all 1200g) : lc"hapin et aLfl2009l) ; ICoppin et al, 
Daddi et al. (2009b a)7 TDevlin et al] 120091 ); IKnudsen et al 



2006) 
200, 
200, 
2009) 
(2010 



selection effect. It is however possible that the correla- 
tion breaks down at even earlier epochs due to changes 
in star formation processes, e.g., suppression of radio 
emission in inver se Compton losses o ff the CMB photons 
as sug geste d by iLacki et all (|2009f ) , lLacki fc Thompson! 
(l200l. andlMurphvl (l2009f)7 

In iMichalowski et al.l (|2010f ) we analyzed the full UV- 
to-radio s pectral energy di s tribut ions (SEDs) of 76 SMGs 
from the Chapma n et al.l (|2005t) sample with spectro- 
scopic redshifts up to z < 3.6. Here we extend that 
study by analyzing the sample of all spectroscopically 
confirmed SMGs at z > 4. The main objective of this 
paper is to characterize the required efficiency of dust 
producers (SNc and AGB stars) at these early epochs 
of the evolution of the Universe. In Section [5] the SMG 
sample is presented. We outline our methodology and 
derive the properties of SMGs in Section [3] and discuss 
the implications in Section [4] Section [5] closes with our 
conclusions. We use a cosmological model with Ho = 70 
km s" 1 Mpc~\ tt A = 0.7 and fl m = 0.3. 

2. SAMPLE 



We selected all 
CO) re dshifts at z 
(2009 1 ICapak et al 
iDaddi et al.l (|2009bl l 
in ' ^ " 



or 



six SMGs with robus t (optical 
> 4 id entified bv ICoppin et ahl 



(120081). ISchinnerer et a" 



and IKnudsen et alj 

,2. 



2008 



2010) 

the ECDF-S (900 arcmin 2 ; lCjmmiLi£aD [2009; 
lOreve et al.l 120091: IWeiletTaTI 120091) . COSMOS (1080 



aremm 



2. 



Scott et al l 120081). GOODS-N (100 



aremm 



2. 



Hugh es et al.l 119981: iBarger et al.l l2000t [Chapm an et al 



2001; Borvs et al.ll2003l : iSerieant et al.ll2003l: iWang et al 









Flux 


Error 




SMG 


z 


(fim) 


(^Jy) 


(MJy) 


Ref. 


LESSJ033229.4-275619 


4.760 


0.352 


0.087 


0.000 


Coppin09 


LESSJ033229.4-275619 


4.760 


0.433 


0.013 


0.000 


Coppin09 


LESSJ033229.4-275619 


4.760 


0.461 


0.039 


0.000 


Coppin09 


LESSJ033229.4-275619 


4.760 


0.538 


0.059 


0.000 


Coppin09 


SMMJ123711. 7+622212 


4.055 


0.365 


0.063 


0.026 


Capak04 



Note. — When the error is equal to zero, the flux column 
denotes the 3a upper limit. Otherwise it denotes formal flux at 
the position of an SMG. The data for SMMJ163555.5+661300 have 
been corrected for lensing magnification of a factor of 5.5. (This 
table is available in its entirety in a machine-readable form in the 
online journal. A portion is shown here for guidance regarding its 
form and content.) 

References, 

120091); 



2004), and Abell 2218 (11.8 arcmin 2 ; IKnudsen et al 
2006) fields. The photometric data are presented in Ta- 
bles [DandEl The data for SMMJ163555. 5+661300 have 
been corrected for lensing magnifica tion of a factor of 
5. 5 obtained bv IKnudsen et al.l (|2010|) using the model 
of lElfasdottir et al.l ()2007D . 

Our sample is not homogeneously selected. Namely, 
some of the sources are bright enough in the optical to al- 
low spectroscopy, whereas redshifts of some of them were 
measured based on CO emission (one even not detected 
at optical wavelengths). In particular, LESSJ033229.4- 
275619 and MMJ100054.48+023435.9 were targeted 
spectroscop ically, because they are U-dropouts sug gest- 
ing z w 5 ([Coppin et al.1 120091 : ICapak et all 120081 ) and 
their optical counterparts (and hence, redshifts) are 
based on detections at radio wavelengths, which likely 
biases them toward high star formation rates (SFRs). 
The CO lines from SM MJ123711.7+622212 (called 
GN20 in lPope et all [20061) and SMMJ123709.5+622206 
(GN20.2a) were detected serendipitously while observ- 
ing an angularly close galaxy at z — 1.522 ([Daddi et al.l 
I2009bl) . and the CO emission of SMMJ123633.8+621408 
(GN10) was searched for under the assumption that 
it is a member of the p rotocluster structur e contain- 
ing GN20 and GN20.2a ([Daddi et al.1 l2009aT ) . Finally, 
SMMJ163555.5+661300 was detected because of its 
lensing magnification and is therefore intrinsically the 
faintest member of our sample. 

It is therefore not easy to quantify how selection effects 
influence our results. Very likely our sample is biased to- 
ward high luminosity objects, i.e., with high SFRs. This 
is supported by the fact that SMMJ163555. 5+661300, 
magnified by lensing and the only member of our sam- 
ple not detected in the radio, has a much lower SFR 
than blank-field members of our sample (Table [3]). How- 
ever, given the significance of the Spitzer IRAC detec- 
tions (rest-frame ~ 1-2 /im) of > 10cx, our sample is not 
biased against low stellar masses. 

3. SED FITTING AND RESULTS 

We applied th e SEP fitting m ethod detailed in 
IMichalowski et al.l (|2008l I2009L I2010L see therein a dis- 
cussion of the derivation of galaxy properties and typical 
un certainties) based on 35 000 templates in the library 
of llglesias-Paramo et al.l (l2007fl. plus some tem plates 
of iSilva et al.l ([1998D an d IMichalowski et all ([2008D . all 
developed in GRASIL ([Silva et al.l I1998D . The tem- 
plates cover a broad range of galaxy properties and 
were test ed to reproduce the SEDs of high-redshift 
galaxies (TSilva et all 119981: llglesias-Paramo et al.l 120071 : 



galaxies (Su va et al.l I199«t llgle 
IMichalowski et al.l 120081 12010). 



Their star formation 
histories are assumed to be a smooth Schmidt-type law 
(SFR proportion al to the gas mass to some power, see 
ISilva et al.l Il998l for details) with a starburst (if any) 
on top of that starting 50 Myr before the time at which 
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10° 10 1 10 2 10 3 10 4 10~ 1 10° 10 1 10 2 10 3 10 4 
Rest Frame Wavelength (jim) 



Fig. 1. — Spectral energy distributions (SEDs) of z > 4 SMGs. Solid lines: the best GRASIL fits. Squares: detections with errors, in most 
cases, smaller than the size of the symbols. Arrows: 3<r upper limits (values marked at the base). The data for SMMJ163555. 5+661300 
have been corrected for lensing magnification of a factor of 5.5. 
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the SED is com puted, mere are seven tree p a 
in the library of llglesias-Paramo et al.l (j2007) : the nor- 
malization of the Schmidt-type law, the timescale of the 
mass infall, the intensity of the starburst, the timescale 
of the molecular cloud destruction, the optical depth of 
molecular clouds, the age of a galaxy and the inclination 
of a disk with respect to the observer. 

We scaled all the SEDs to match the data (detections) 
and chose the one with the lowest \ 2 ■ In an but one case 
(see below) the upper limits did not provide a constraint 
on the best-fitting model. 

In all b ut one case we obtained the b est fits using the 
library of llglesias-Paramo et al.l (|2007j ). For GN10 we 
fitted a model correspo nding to a 0.1 Gy r old progenitor 
of an elliptical galaxy (Sil va et al.lll998l) with modified 
maximum grain temperatures from 400 K to 100 K in 
order to suppress strong mid-IR emission in the original 
model, otherwise inconsistent with the data. All other 
models failed to reprod uce its extremely red observed- 
frame 2.1-3.6 ^irn color (jWang et al.ll2009f ). 

The best fits 3 are shown in Figure [T] It is apparent 
that GRASIL models can reproduce the SEDs of even 
such distant galaxies. The resulting properties of the 
galaxies are listed in Table 

In order to assess the uncertainties of the derived pa- 
rameters, we calculated the confidence intervals using 
SED models with A\ 2 < 1. The resulting uncertain- 
ties are a factor of ~ 2-4 for SFRs and Ljr; a factor of 
~ 1.5-2.5 for M*; a factor of ~ 2-5 for M burst ; - 10 K 
for Td; and 0.3-1.0 mag for Ay- The uncertainty on Td 
introduces a factor of ~ 1.5-2.5 error on Md and allow- 
ing the emissivity j3 to vary in the range 1-2 results in 
an uncer t ainty of a factor of < 3 on Md- The choice of a 
iSalpeter] (|1955l ) initial mass function (IMF) with cutoffs 
of 0.15 and 120 M Q introduces a maximum systematic er- 
ror of a factor of ~ 2 in the deter mination of the stellar 
masses and SFRs (|Erb et al.ll200H . 
Our estimates of SFRs, Ljr and Md are con- 
within a factor of < 3 wit h those ob- 

by |Corjrjin et al.l (120091) . iCapak et all 
Schinnerer et al.l (I2008D. lYounger et al.l 



sistcnt 
taincd 



ICasev et al.1 



Dad di et all (|2009bl laf) . 



iKnudsen et alJ (|2010T >. and lWang et al.1 (|2009f ) (compare 
Tables [3] and IH. Ou r M* es timates are consistent with 



those of lDaddi et ail (|2009bl lah. and IWang et all (120091) 
after t aking int o acco unt that iDaddi et al.1 ()2009bl ial) 
used a IChabrieil ((2003) IMF res ulting in st ellar masses 
1.8 times lowe r than for the ISalpeterl (|1955f ) IMF 
(|Erb et al.M2006h . 

However we obtained stellar masses for 
LESSJ033229.4-275619 and M MJ100054.48+ 23435 9 

~ 10 times la rger than ICoppin et al.1 (j2009) 4 

and ICapak et al.1 ([2008), respectively, and for 
SMMJ163555.5+66130 - 5 times larger than 
IKnudsen et~all (|2010t l. The differ ence can be ex- 
plained by the fact that ICoppin et al.l (|2009l ) assumed a 
mass-to- light ratio of MJLk = 0.1 and the determina - 
tions of ICapak et all (|2008( 1 and IKnudsen et"all (|2010f ) 
correspond to M*/Lk ~ 0.03 and ~ 0.17, respectively. 
These are very low values compared to lower- 

3 The SED fits can be downloaded from 
http : / /archive . dark- cosmology . dk 

4 Note that our estimate agrees with that of Stark et al. (2007). 
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TABLE 4 

Properties of z > 4 SMGs Derived by Other Authors 





SFR 


Lm 


Af» 


SMG 


(M@ yr" 1 ) 


(lO^L©) 


(10 11 M e 


LESS 


1000 


0.61 


0.5 


MM 


1000-4000 


0.5-2.0 


0.1 


GN20 




1.0-5.0 


2.3 


GN20.2a 




0.5-1.6 


0.5 


GN10 


2400 


1.2-2.5 


1.0-3.0 


SMM 


230 


0.13 


0.16 



0.5 



Ref. 

Coppi n et al. (20091 

Capak et al. (2008) . Schinnerer et al. (2008) 

Daddi et al. (2009b) . Younger et al. (2008b) , Casey et al. (2009) 

Daddi et al. (2009b) . Casey et al. (2009) 

Daddi et al. (2009a) . Wang et al. (2009) 

Knudsen et al. (2010) 



Note. — The first column lists all SMGs in the order given in Table [3] Only the first parts of their names are given for brevity. 



redshif t galaxies (e.g . . iDrory et al.ll2004l; iPortinari et all 
20041: lLabbe et alf 120051: ICastro Ceron et all l200€' 



20M Ivan der Wei et all 12005 iCourtv et all 120071 : 
Michalowski et all |2010T ISavaglio et all 120091) giving 
lower limits on stellar masses. On the other hand, we do 
not assume mass-to-light ratios, but derive them from 
the stellar population models incorporated in GRASIL. 



4.1. 



4. DISCUSSION 
Formation of Stars in z > 4 SMGs 



All z > 4 SMGs in our sample are characterized 
by an extremely strong starburst episode (average SFR 
~ 2500 Mq yr" 1 , Column 5 of Table EJ) during which 
a substantial fraction (average ~ 45%, Column 9) of 
their stellar population was formed. They are there- 
fore manifestations of the strongest known star-forming 
events in the universe. Their high stellar mas ses (average 
~ 3.6 x 10 11 M Q ) agree with a suggestion of iDave et all 
(2009) based on numerical simulations that the most 
rapidly star forming galaxies coincide with the most mas- 
sive galaxies. 

The pr operties of the galaxies are within the ranges 
found bv lMichalowski et all (120101) f or the 1.7 < z < 3.6 
SMG sample (jChapman et alT T2005). However, z > 4 
SMGs tend to be among the most star-forming, least 
massive, and hottest SMGs and exhibit the highest frac- 
tion of stellar mass formed in the ongoing starburst. 
Namely, 43% of 1.7 < z < 3.6 SMGs have lower 
SFRs than any of the z > 4 SMGs 5 , whereas only 2% 
have higher SFRs; 30% of 1.7 < z < 3.6 SMGs have 
higher stellar masses, whereas only 8% have lower stellar 
masses; 55% have a lower fraction of stellar mass formed 
in the ongoing starburst episode, whereas only 2% have 
higher fraction; and 28% have lower dust temperatures, 
whereas only 4% have higher temperature. This can 
be interpreted as SMGs at z > 4 representing earlier 
stages of the formation of submillimeter-bright galaxies 
at which the pre-existing stellar population is less pro- 
nounced and therefore the ongoing starburst episode con- 
tributes more to the final stellar mass. 

However, we note that we cannot exclude the possibil- 
ity that these galaxies are AGN dominated, suggested 
by their high IR luminosities and dust temperatures. In 
that case the SFRs and M„ we derive would be upper 
limits. 

If GN10 is not AGN dominated, then the fact that 
we could only fit its SED using a template correspond- 
ing to a young (0.1 Gyr old) progenitor of an elliptical 

5 We exclude SMMJ163555.5+661300 from this analysis, be- 
cause its unlensed submillimeter flux of ~ 2 mjy makes it impos- 
sible to be det ected by SCUBA in th e blank-field survey similar to 
those used by Chapman ct al. (2005). 



galaxy supports the hypothesis that SMGs evolve into 
ellipticals. According to our model, after 1.5 Gyr GN10 
will evolve into a massive elliptical containing ~ 10 12 Mq 
stars. In order to explain its extremely red observed- 
frame 2.1-3.6 /jm color, we did not need to invok e an ol d 
stellar population as sugg ested bv IWang etHTl (2009). 
Its spectrum is reproduced by a young stellar population 
residing in molecular clouds heavily obscured by dust 
with Ay ~ 7.8 mag (Tabled Column 14). This esti- 
mate is similar to t he values obtained u sing only near-IR 
and op tical data bv lDaddi et all (|2009af ) and lWang et all 

On the other hand the remaining z > 4 SMGs are 
only moderately obscured (Ay ~ 2 mag; Column 14). 
We checked that this value of extinction is roughly con- 
sistent with the expectation based on the estimation of 
the amount of energy reprocessed by dust. A crude 
estimate of the UV (0.2 /mi) extinction can be calcu- 
lated using Ao^m = 2.5 log(SFRi R /SFR uv ). This gives 
Ao.2fj,m. ~ 20 mag for GN10 and ~ 3.8-6.3 mag for the 



other objects . Ass uming that A0.2 



3 x A v (e.g. 



iGordon et all [2003). this corresponds to Ay ~ 6.6 mag 
for GN10 and Ay ~ 1.2-2.1 mag for the other objects. 
These values are consistent with those derived from the 
SED modeling (Table [3]). 

We estimated the SMG comoving volume densities as 
a sum of SFRs (or Lir, or stellar or dust masses) of z > 4 
SMGs 6 divided by the comoving volume in the redshift 
interval 4-5 (5.9 x 10 6 Mpc 3 ). We obtained p S FR = 
2.6 x 10~ 3 Mf7) yr" 1 Mp c~ 3 (a similar value was found 
by lYounger et all l2009al based on candidates for high- 
redshift SMGs), logp L iR = 7.17L© Mpc~ 3 , logp* = 
5.54 Af Mpc -3 , and logpdust = 2.81 Mq Mpc" 3 . How- 
ever one must keep in mind that these numbers could be 
affected by cosmic variance, because three out of five 
SMGs co nsidered here are m embers of a protocluster 
structure (Da ddi et al.l 2009b. al The contributions to 
the cosmic SFR and stellar mass densities of SMGs at 
these redshifts are 4.4% and 1.0%, respectively, using 
the co mpilation of the total values in Michalowski ct al. 
(|2010l Tables A. 4 and A. 5). These numbers indicate that 
currently detected SMGs did not contribute significantly 
to the cosmic star formation history at z > 4, but our 
estimates should be regarded as lower limits since more 
of such distant SMGs could still be undetected in the 
fields discussed here. 

4.2. Producers of Dust in z > 4 SMGs 

The dust masses we find for z > 4 SMGs (Ta- 
ble [3l Column 11) are similar to those derived for 

6 We again exclude SMMJ163555. 5+661300 because it is lensed. 
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TABLE 5 

Dust Yields Per Star Required to Explain Dust in z > 4 SMGs 











Dust Yields (M Per St 


ar) 




Dust Producer 


IMF 


Total Mass 


LESS 


MM 


GN20 


GN20.2a 


GN10 


SMM 


AGB (2.5 < M < 8M ) 


Salpctcr 




0.03 


0.06 


0.04 


0.17 


0.05 


0.07 


AGB (2.5 < M < 8M e ) 


Salpeter 


M, - A/burst 


0.04 


0.31 


0.06 


0.20 


oo 


0.09 


SN (M > 8M @ ) 


Salpeter 




0.14 


0.22 


0.15 


0.65 


0.18 


0.27 


SN (M > 8M ) 


Top-heavy 


M* 


0.05 


0.08 


0.05 


0.23 


0.06 


0.10 



Note. — The IMF is either Salpctcr (1955) with a = 2.35 or top-heavy with a = 1.5. The total mass indicates if the entire stellar mass 
(M») was used to calculate the number of stars (Equation lAll l. or if stars created during the ongoing starburst were excluded (M* — A/burst)- 
The last six columns contain dust yields for all SMGs in the order given in Table|3] Only the first parts of their names are given for brevity. 



4-6 quasars (a fewxlC^Mp 



Dunlop et all 119941: 
Benford et al.lll999t [Archibald et alj|2001t lOmont et al 
20011: iPriddev fc Mcrvj a horj 120011: iPriddey et al 
20031 120081: Jlsaak et al.1 l2002t iBertoldi et al.1 12003 



Robson et alj|2004: iBeelen et alj |2006l: IWang : et al.ll2008t 
Martmez-Sansigre et alj|2009|) . These huge dust masses 
indicate that dust in SMGs was efficiently formed and 
able to survive even when the age of the universe was 
only 1.2-1.5 Gyr. 

As detailed below we find that i) AGB 
stars are efficient and fast enough to form 
dust in the LESSJ033229.4-275619, GN20 and 
SMMJ163555. 5+661300; it) only SNe are efficient and 
fast enough to form dust in MMJ100054.48+023435.9, 
GN20.2a, and GN10, as long as the unusually high dust 
yields derived for Cassiopeia A and Kepler SN remnants 
are correct and typical, or if a top-heavy IMF and/or 
significant dust growth in the interstellar medium (ISM) 
is assumed. If these assumptions are correct, then SNe 
could also produce dust in SMGs mentioned in i. 

These conclusions are not significantly affected by the 
uncertainties reported in Section [3] Even if dust masses 
were higher or stellar masses lower by a factor of a 
few, the required dust yields per SN woul d not exceed 
that of Cassiopeia A ([Dunne et al.l l2009f ) and Kepler 
(jGomez et al.| [2009). It would however be more critical 
for AGB stars if dust masses were higher than we de- 
rived, since the required dust yields for LESSJ033229.4- 
275619, GN20, and SMMJ163555. 5+661300 are already 
on the high end of the theoretically allowed yields. How- 
ever, a decrease of dust masses by a factor of a few would 
not change our conclusion that AGB stars are not effec- 
tive enough to form dust in MMJ100054.48+023435.9, 
GN20.2a, and GN10. 

4.2.1. Asymptotic giant branch stars 

In order to investigate whether AGB stars can be re- 
sponsible for dust production in z > 4 SMGs, we es- 
timated (see the Appendix) the average dust yields re- 
quired per star with mass 2.5 < M < 8 M@ and main- 
sequence lifetime in a range 1 Gyr-55 Myr (calculated 
as 10 10 yr x [M/M p 2 ' 5 ; e.g.. iKippenhahn fc Weigertl 



ll99Cf l. The lower mass limit was chosen to ensure that 
the stars considered can start producing dust within the 
age of the universe at the redshifts of our sources (1.2-1.5 
Gyr). 

The results are listed in the first row of Table [5] They 
are independent of the assumed star formation history 
of galaxies, but depend only on derived dust and stel- 
lar masses, assumed IMF and measured redshifts. We 
find that each AGB star would need to produce ~ 0.03- 
0.07 M Q of dust in order to explain the dust in z > 4 



SMGs, excluding GN20.2a (see below). These num- 
bers are cl ose to the highest theoret i cal dust yields of 
AGB stars ([Morgan fc Edmundjl200llFerrarotti fc Gaill 
2006, the total dust yields per AGB star are approxi- 
mately metallicity independent) making them plausible 
dust producers. Assuming a top-heavy IMF does not 
change this result significantly (the required dust yields 
increase by ~ 30%), because the choice of the slope of 
the IMF affects the number of high- and low-mass stars, 
leaving the number of intermediate-mass stars approxi- 
mately constant. For GN20.2a the required dust yield 
(0.17 Mq) is too high to claim that AGB stars formed its 
dust. 

However, MMJ100054.48+023435.9 and GN10 formed 
the majority of their stars in the ongoing starburst 
episode (Column 9, Table [3]), which is too short for the 
2.5 < M < 8M Q stars considered above to finish their 
main-sequence phase. Therefore AGB stars could not 
contribute to the dust production in these two galax- 
ies. To quantify this, we calculated the required dust 
yields for AGB stars taking into account only stars that 
were born before the ongoing starburst (replacing M* 
by M* — A/burst in Equation IA1[) . The resulting yields 
(the second row of Table© for LESSJ033229.4-275619, 
GN20, GN20.2a, and SMMJ163555.5+661300 do not dif- 
fer significantly from our previous estimate (the first row 
of Table [5]), because in this way we removed < 30% of 
stars (those formed during the ongoing starburst). How- 
ever the yields for MMJ100054.48+023435.9 and GN10 
become too high to claim that AGB stars formed dust in 
these galaxies. 

A potential limitation of this claim is the uncertainty 
in determining the fraction of stellar mass formed during 
the ongoing starburst episode. If more stars were formed 
in the past, then AGB stars could be responsible for dust 
production in these two galaxies. However, even if it 
was the case, then the current SFRs would be the same as 
we derive, because they are fixed by strong submillimeter 
emission. Then the ongoing starburst becomes unfeasi- 
bly short (< 25 Myr) in order not to produce more stars 
than is inferred from the optical to near-infared part of 
the spectra. Therefore, more stars in these two SMGs 
could have been formed before the ongoing starburst only 
if the current SFRs are overestimated due to a significant 
AGN contribution. 

4.2.2. Supernovae 

We repeated the analysis for SNe, i.e., calculated the 
required dust yields per one massive star (> 8M@). 
The yields are ~ 0.15-0.65 M of dust per SN (the 
third row of Table [5]), consistent with the theoretical 
works (though without dust grain destruction implc- 
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2.5 



2.0 -, 



1 .5 - 




• 4.0<z (This work) 
2.6<z<3.6 (M09) 
2.1<z<2.6 (M09) 
1.4<z<2.1 (M09) 

• 0.1<z<1.4 (M09) 
□ averages 

Local Value 

Scatter 

Linear fit 







2 3 
Redshift 



mente d) of iTodini fc Ferraral (|2001h and iNozawa et all 
(|2003l) ; with a value predicted by iDwek et al.l (|2007| ) 
to account for dust in a z ~ 6.4 quasar; and with 
su bmillimeter observation o f SN remna nts Cassiopeia 
A (iDunne et al.ll2003l [20 09) and Kepler (|Morgan et al.l 
I20031 iGomez et al.l l2009f K There is a debate about 
the la t ter r es ults on Cass i opeia A and Kepler (e.g. , 
Dwell [2001 iKrause et al.l [200l IGomez et all [200a 
Wilson fc Batrlall2005HBlair et alj|2007t ISibthorpe et al l 



2009J), but if thedust yields of IDunne et all (|2009t) and 



Fig. 2. — Ratio of the infrared (8-1000 /im) and radio luminosities q = log(Lia/3.75 X 10 12 / I1.4) as a function of reds hift of SMGs. Both 
SMGs at z > 4 (red circles) and at lower redshifts (yellow, green, blue, violet circles, divided into four redshift bins) from Michalows ki et all 
(2010, M09) are shown. In the redshift range 1.4 < z < 5.0 no significant evolution of the IR-radio correlation is found for SMGs. The 
average q for SMGs (squares) is howe ver offset toward higher radio luminosities (factor of ~ 2.1-2.3) from the local relation (q = 2.64: solid 
line with scatter, 0.26: dotted lines: lBeTill2003I V A linear fit to the data (dashed line) resulted in q decreasing with redshift, but only at 
~ 2cr significance. 

have b een claimed to reproduce the number counts of 
SMGs: lBaugh et alJl200alFontanot et al.ll200l . Chang- 
ing the IMF slope from a = 2.35 to a = 1.5. consi stent 
with values for low-mass star clusters (IScalol I1998T) an d 
a limit derived for a proto-star cluster (|Sternberdll998|l 
resulted in the required dust yield decreasing to 0.05- 
0.23 M Q (the fourth row of Table EJ. 

The second possibility is that SNe provided only the 
dust seeds and that the bulk of the dust mass w as accu- 
mula ted during grain growth in the ISM (e.g., iDraind 
l2003h . The timescale of this process is typically less 
than a few x 10 M yr (|Hirashital 120001: IZhukovska et all 
120081 : lDraindl2009l ). i.e., short enough to contribute sig- 
nificantly to the growth of dust mass in z > 4 SMGs. 
However, the extreme environments of z > 4 SMGs may 
make it difficult for grains to grow before they are de- 
stroyed. The hypothesis of dust growth can be tested 
by investigation whether the grains formed by SN rem- 
nants are smaller than those present in z > 4 SMGs. 

4.3. The IR-Radio Correlation at z > 4 

Despite the differences compared to the 1.7 < z < 3.6 
population fSection l4.1[) . the mean IR-to-radio luminos- 
ity ratio (q = log[LiR,/3.75 x 10 12 /ii. 4 ] with radio K- 
corrections assuming a slope of —0.75) for z > 4 SMGs 
of 2.32 ± 0.20 is con sistent with the values d erived for 
lower-redshift SMGs (|Michalowski et alll2010t ). At z > 4 
we find an offset of Aq ~ —0.32 (a factor of ~ 2.1) from 



Gomez et all ((2009!) were correct and typical, then SNe 
would be efficient enough to account for the dust in all 
z > 4 SMGs. 

However our, required estimates are at least an 
order of magnitude higher than any other observed SN 
dust yields , which are t ypical ly in the range ~ 10~ 3 - 
10~ 2 M^ (TGreen et all l200l fe orkow ski et all 120061: 



Sugerman et alj|200a lErcolano et alj|2007t iMeikle et al 



20071 iRhoetall I2008L 120091: 



Kotak et al.| __|2009: 



Sakon et all 120091: ISandstrom et al.H2009t IWesson et all 



2009( 1 and theoretically pr edicted dust masses able t o 



survive in SN remnants dBianchi fc Schneider! 120071 ). 
This apparent difficulty in explaining dust production 
in z > 4 SMGs can be resolved with a combination of 
two plausible effects. 

Approximately half of the discrepancy can be ac- 
counted for with a top-heavy IMF giving more SNe per 
unit stellar mass (both top-heavy and Salpeter IMFs 
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the local value of q = 2.64 (|Bellll2003l ). though due to 
the small sample, this offset is significant only at ~ 1.6cr 
level. 

Hence, wit h our mul t i- wav elength approach we confirm 
the results of lMurp hv (2009), who derived the mean q = 
2.16 ± 0.28 for z > 4 SMGs. Moreover, our estimate is 
consistent with a value of 2.4 1 ± 0.20 derived for 250 fxm- 
selected galaxies at z < 3 bv llvison et al.l (|2010t ). 

The offset of Aq ~ —0 .32 i s consistent with a hypoth- 
esis of lLacki et aTJ(|2009| ) and lLacki fc Thompson! (|2009h 
that SMGs are "puffy starbursts" (vertically and radially 
extended galaxies with vertical scale heights ~ 1 kpc) ex- 
periencing weaker bremsstrahlung and ionization losses 
resulting in stronger radio emission. We cannot however 
exclude an AGN contribution boosting their radio fluxes. 
Since the redshifts of four out of six z > 4 SMGs have 
been measured independently of radio detections, the ra- 
dio excess cannot be a result of a bias against radio-faint 
sources. 

In Figure [2l we show the q- values as a function of red- 
shift for both z > 4 SMGs discussed in this paper (red 
circles ) and lower-redshift SMGs from Michalowski et al. 
(|201Ct ). The IR-radio correlation of SMGs does not 
show any significant evolution in the redshift range 
1.4 < z < 5.0. A linear fit to all the data results 
in a low-significance (~ 2a) ind ication that q dec reases 
with redshift (as also found by llvison et al.l 120101 from 
a stacking analysis of 24 /im-selected galaxies) in the 
form q = (-0.078 ± 0.038)z + (2.50 ± 0.09). To date, 
there are only two q determinations at higher redshifts 
than presented here, namely z = 6.2 and 6.42 quasars 
(|Carilli et al.ll200l iBeelen et al.ll200l q = 1.8-2.2). 

4.4. Gas-to-Dust Ratio at z > 4 

Using the molecular gas mass e stimates (based 
on C O [4-3] line observ at ions) from ISchinnerer et al.l 
((20081) and iDaddi et al.l (l2009bl lah we derive gas-to- 
dust ratios of M H2 /M d = 73, 47, 22, and 98 for 
MMJ100054.48+023435.9, GN20, GN20.2a, and GN10, 
respectively. We note that these values are affected 
by significant uncertainties of the conversion from CO 
line strength to H2 mass. The mean value of 60 is con- 
sistent with 54+!i estimated for z ~ 1-3.5 SMGs by 
iKovacs et al.l (|2006l ) using the CO survey of iGreve et al.l 
(12005ft . SMGs at z > 4 have one of the lowest gas-to- 
dust ratios compared to galaxies selected at other wave- 
lengt hs, e.g., the Milky Way (~ 90 - 400;ISodroski et al. 
19971) other spirals (~ 1000 ± 500; iDevereux fc Young 
19901 : iStevens et al.ll2005l ) , the nuclear regions of local lu- 
minous IR galaxies , ultraluminous IR galaxies (120 ±28; 
IWilson et al.ll2008D. and of local, far-IR-selected galaxies 
(~ 50: ISeaquist et al.ll2004l ) (all these results are based on 
IR and CO data, and are therefore directly comparable 
to our estimates). This is not surprising, since SMGs are 



selected by submillimeter emis sion. Similarly, a low y alue 
of Mh 2 /Mcl — 30 was found bvlBertoldi et all (l2003bh for 
z = 6.42 quasar, but ICox et alT i 20021 ) reported a higher 
value of Mh 2 /M,i = 150 for z ~ 4 quasar. The small 
number of SMGs and quasars with derived gas-to-dust 
ratios hampers a comparison of these two samples. 

5. CONCLUSIONS 

We have analyzed the SEDs of six spectroscopically 
confirmed z > 4 SMGs. Our results provide constraints 
on dust production at these early epochs of the evolution 
of the universe. We find that AGB stars are efficient and 
fast enough to form the dust residing in three of these 
galaxies. However, for the three remaining SMGs only 
SNe are efficient and fast enough. The high required SN 
dust yields hint at a possibility that their stars may be 
distributed according to a top-heavy IMF and/or that 
dust grains are substantially grown in the ISM. Since 
the majority of the stars in these galaxies were formed 
on very short timescales, these properties are very likely 
to be similar to those of the first galaxies beyond redshift 
6, which had been building up their stellar populations 
only for several hundred Myr, which elapsed since the 
big bang. 

We present evidence that the IR-to-radio luminosity 
ratios of z > 4 SMGs are consistent with that of lower- 
redshift SMGs and are offset from the local relation by 
a factor of ~ 2.1. 

A comparison of the z > 4 SMGs with the lower- 
redshift sample, in particular their high SFRs, dust tem- 
peratures and fraction of stars formed during the ongoing 
starburst as well as low stellar masses, reveals that we 
start to see SMGs at earlier stages of their evolution. 

The improved mapping speed and sensitivity of the 
new SCUBA2 camera will enable studies of the evolu- 
tionary sequence of SMGs using much bigger and more 
homogeneously selected samples. Moreover, the study of 
dust production presented here will be pushed forward 
with a synergy of Herschel, SCUBA2, and ultimately 
ALMA. These facilities will provide a broad wavelength 
coverage at the IR, which will allow accurate determina- 
tion of dust temperatures and, in turns, its mass. 
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viding his SED templates. 
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APPENDIX 



DUST YIELD CALCULATION 



We calculated the dust yield per star required to explai n dust mass in a galaxy in the following way. In an IMF 
with Mmin = 0.15, A/max = 120 Mq, and a slope a = 2.35 (|Salpeterlll955l or a — 1.5 for top-heavy IMF), the number 
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of stars with masses between Mo and Mi in the stellar population with a total mass of M* can be expressed as 

/■Mi 

/ M- a dM 

N(M <M <M 1 )=M*-^L. , (Al) 

/ " M- a MdM 

where the denominator provides a normalization so that a total mass is equal to M* . For SNe we assumed Mq = 8 M© 
and Mi = M max = 120 M , whereas for AGB stars: M = 2.5M and Mi = 8M Q . 
The average dust yield per star is equal to the dust mass divided by the number of stars, N(M < M < Mi). 
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